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I .  Summary 

During  the  first  year  of  support  of  AFOSR  Grant  85-0253,  Dr.  G.  Chen  and 

his  collaborators  made  continuous  progress  on  the  resear  of  control  and 
stabilization  problems  in  structural  dynamics.  His  primary  emphases  were  on 
the  modelling,  designs,  placement  and  analysis  of  stabilizers  and 
controllers.  Two  primary  types  of  vibrating  structures  were  treated:  the 
second  order  wave  equation  such  as  strings  and  cables,  and  the  fourth  order 
equation  such  as  beams.  For  those  structures,  he  was  able  to  classify  the 
types  of  stablizers,  analyze  their  damping  behavior,  and  determine  the 
mechanical  designs.  Various  methods  such  as  characteristics,  asymptotic 
estimation,  energy,  functional  analysis  and  the  Legendre  spectral  method  were 
used  to  study  the  problems. 

Some  relevant  theory  and  methods  for  general  distributed  systems  and 
equilibrium  problems  were  also  developed. 


II .  Research  Program 

The  research  objectives  of  our  research  program  last  year  consisted  of 
the  study  of  the  following 

i)  Asymptotic  spectral  distribution  and  practical  implementations  of 
various  damping  conditions  at  boundary  or  joints  for  coupled  strings  or 
beams : 

ii)  Numerical  computation  of  the  spectrum  of  a  boundary  damped  dynamic  bi- 
harmonic  equation  modelling  a  plate: 

iii)  Stabilization  and  control  of  nonlinear  vibrating  systems. 

So  far,  we  have  already  completed  the  first  phase  of  these  objectives,  as  can 
be  seen  from  our  list  of  publications. 

During  the  second  year  of  support,  we  plan  to  enter  the  second  phase  of 
our  research  program  by  continuing  the  investigation  of  the  above,  focusing  on 
more  complicated  systems.  Meanwhile,  we  are  also  planning  to  research  the 


following  new  objectives: 


:'vdes 

j  /v... .  j,  o  /  or 

I  Id) 


Ci’py 
.  »Nsr*  r  * i 


/9y 


a 


ft 


/  .*  .  r%  r  ,  *  .  r«  J*  ‘  V’  .  ’  .  '  '  •  *  •  *  *."  <*  *  ' 


-2- 


iv)  Analysis  and  control  of  dynamic  frame  structures; 

v)  Study  of  new  beam  models  such  as  the  one  recently  proposed  by  D.  L. 
Russell  which  exhibits  better  damping  characteristics. 

vi)  Higher  dimensional  membrane  and  plate  equations  with  nonlinearity. 
They  will  constitute  the  major  portion  of  our  research  program  for  the 
second-year. 

III.  Research  Accomplishments  and  Work  in  Progress 

We  first  describe  our  research  findings  under  headings  [ A ] - [D ]  below: 
[A]  Analysis  and  designs  of  stabilizers  for  second  order  systems-vibrat ing 

strings  and  cables . 

The  standard  model  is  the  wave  equation 

m  -  t  =  0>  0  <  x  <  L.  t  >  0. 

at*2  ax'2 

where  y(x,t)  represents  the  vertical  displacement  at  position  x  at  time 
t,  m  is  the  mass  density  per  unit  length,  and  T  is  the  tension 
coefficient.  The  energy  of  vibration  at  time  t  is 

E(t)  =  1  JQ  [■y^ix.t)  +  Ty  \  ( x ,  t )  ]  dx . 


If  a  stabilizer  is  placed  at  the  boundary  x  =  L,  then 
(1)  ^  E(t)  =  Tyt(L,t)yx(L,t)  <  0. 

assuming  that  the  boundary  condition  at  x  =  0  is  energy  conserving.  For 
(1)  to  hold,  we  must  use  proportional  control 


» 
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(2)  Ty x ( L , t )  =  -k2y t ( L . t ) .  k2=feedback  gain  >0,  for  all  t  >  0, 

i .  e .  , 

force  is  negatively  proportional  to  velocity  at  x  =  L. 

Therefore  a  standard  viscous  damper  is  sufficient  for  this  purpose,  as  shown 
below: 


F  i  g,  1  A  vibrating  string  with  a  damping  device  at _ le_f_t — end 

The  above  fact  is  well  known.  Furthermore,  if  one  chooses  the  feedback  gain 
to  be 


k2  =  T  v/T7m  . 

then  the  damper  action  (2)  becomes  a  so  called  characteri Stic  impedance 
condition,  which  causes  maximum  energy  loss  to  the  vibrating  system. 

Now,  let  us  consider  the  case  when  the  stabilizer  is  installed  at  an 
in-span  point,  say  at  x=a.o<a<L.  Assume  that  the  boundary  conditions 
at  x  =  0  and  x  =  L  are  energy-conserving.  Then 


(3)  [Total  energy  at  time  t]  =  (■§  J  [my^  ( x .  t ) -*-Ty2  ( x ,  t )  ]dx 


\  J  [my2 (x . t )  *  Ty2 ( x , t ) ]dx }  =  Tyx(«' . t )yt («" , t ! 


Tyx(<*  ,t)yt(a  ,tj  <  0. 


These  are  only  two  possibilities  which  can  cause  the  negativity  of  the  right 
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hand  side  of  (3)  : 

(Al)  Continuity  of  displacement  (and  velocity) 
y(«",t)  =  y(«\t) 
discontinuity  of  vertical  force 

Tyx(a‘,t)  -  Tyx(a*.t)  =  -k jyt (a- , t )  ,  kj=feedback  gain  >  0; 

(A2)  Discontinuity  of  displacement  and  velocity 

yt(«~.t)  -  yt(a’.t)  =  k2Tyx(a~ '  t)  >  k|=feedback  gain  >  0, 

continuity  of  vertical  force 

Tyx(«-.t)  =  ryx(«-.t). 

The  mechanical  designs  of  stabilizers  of  the  above  two  types  are  given 
in  Figures  2-5. 

The  mathematical  analysis  of  damping  behavior  associated  with  the  above 
designs  is  given  in  our  paper  [2],  along  with  numerical  results. 

[B]  Analysis,  designs  and  placements  of  stabilizers  for  fourth  order 
systems-beams 

These  include  the  Euler-Bernoul 1 i ,  Rayleigh  and  Timoshenko  beam  models. 
We  will  use  the  Euler-Bernoull i  beam 

m  *2y 4*11  +  EI  _  o,  0  <  x  <  L.  t  >  0. 

at  uxJ 

to  illustrate  our  designs  and  analysis.  In  the  above,  y(x.t)  is  the 
vertical  displacement,  m  is  the  mass  density  per  unit  length,  and  EI  is 
the  flexural  rigidity  coefficient.  The  energy  of  vibration  at  time  t  is 

E  ( t )  =  \  j^[myj(x.t)  *  dy^x.tlldx. 


♦ 


t  S 


t 

i 

Again,  let  us  install  a  stabilizer  at  the  right  end  x  =  L,  and  assume  that  ! 

the  boundary  conditions  at  x  =  0  are  energy  conserving.  Then 

(4)  af  E(t)  =  EI[yxx(L.t)yxt(L,t)  -  yxxx( L , t )yt (L , t ) ]  <  0. 

Now,  to  achieve  the  negativity  of  the  right  hand  side  above,  there  are  many 
different  boundary  conditions  and  corresponding  designs  which  can  serve  this 
purpose.  The  following  is  a  more  or  less  exhaustive  list  of  them: 


-El  yxxx(Llt)  =  kj  >  0 

-EIyXx< L ■ t)  =  0 


EIyxxx*E,t^ 

-EIyxx(L,t)  =  k|yxt(L.t).  k£  >  0 


(7) 


yx(L.t)  =  o 

'EIyxxx(Lit)  =  _kiyt(L, t* •  ki  >  0 


(8) 


'y(L.t)  =  o 

-EIyxx(L,t)  =  k|yxt(L.t).  k\  >  0. 


f_E  1  yVvv  ( L  ■  t )  =  -k?yt  ( L .  t )  c.y  (L.t)  .  k2  >  0 

l'EIyxx(L’t)  =  k2yxt(L,t)  +  c2yt(L.t).  k\  >  0. 


where  in  (9).  Cj  and  c2  are  real  constants  satisfying 


» 


i 
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(C1-c2 )f'*J  -  -  k|n2  £  0.  for  ail  <\n  6  R. 

Note  that  throughout  < 5 ) - < 9 ) .  all  terms  have  physical  meanings: 


yt  =  velocity 

yx  =  rotation,  yxt  =  angular  velocity, 

-Elyxx  =  bending  moment, 

-EIyXxx  -  shear 

Thus  the  feedback  schemes  ( 5 ) - ( 9 )  are  all  proportional  controls.  Their 
mechanical  designs  are  indicated  in  Figures  6-10,  cf.  [3]. 

The  analysis  of  eigenf requencies  of  vibration  has  also  been  given  in 

t3j. 

[C]  Stability  criterion  for  point  stabilizers  on  coupled  quasilinear 
vibra-ing  strings 

Consider  two  quasilinear  wave  equations 


a2 

2—5  w(x , t ) 
3t^ 


=  0, 


0  <  x  <  1 , 


(10) 


t 


^—5  w(x. t ) 
dt* 


3x 


0, 


1  <  x  <  2 . 


where  Oj  and  a2  are  nonlinear  functions  satisfying  o^(0)  =  0,  <7  •  ( 11 )  >  0 

for  i  =  1,2.  The  boundary  condition  at  x  =  0  is 

(11)  0 1 ( wx ( 0 • t ) >  -  k2wt(0,t)  =0,  k2  >  0. 


and  at  x  =  2  is 


(12) 


w(2,t) 


0. 
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The  two  strings  are  coupled  at  x  =  1  through 


(13) 


o1(wxd  ,t))  =  o2(wx(  l"  ,  t ) ) 
wt(l'.t)  -  wt(l  +  .t)  =  -kjcr1(wx(r.t)) 


kl  >  o. 


Using  some  fairly  recent  theorems  on  quasilinear  hyperbolic  PDEs,  we  are 
able  to  show  that  if  kg  and  kj  in  (11)  and  (13)  are  nonnegative,  then  the 

solution  decays  exponentially  in  the  C1-norm: 


-at 


(14)  II  w  ( -  .  t )  ||  .  +  II  w(  •  t  )||  .  <  Ke  “ L ,  K,  a  >  0  for  all  t  >  0, 

c‘(0,l)  C"1  (1.2) 

provided  that 


1  ♦  *  «2  *  KjKgK^  -  4CjC2A  "Kj  >  0 


2(1-K1K2K3  ♦  4CjC2A  wK3)  >  0 


1  -  KjK3  -  K,  <■  KjK2K3  -  4  c^A  ^K3  >  0 


are  satisfied,  where 


ci  9  yu^ToT  ,  i  =  1,2.  Cj  > 


A  *  Cj  -  c2  -  k^ctc2 


K|  *  A  (c^  c'>^^icjco) 


K2  ”  A  I c 1  c2  'k 1 C 1 c2 I 


k3  =  I ci  k0i'  (crko> 


for  solutions  of  ( 10)— ( 13)  with  sufficiently  smooth  small  initial  data.  Note 


« 
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that  if  =  0  in  (11),  then  generally  (14)  does  not  hold  no  matter  how 

>  0  is  chosen.  This  shows  that  installing  a  stabilizer  in  the  middle  of 
the  span  (without  installing  another  one  at  the  boundary)  is  not  robust  with 
respect  to  stability. 


[D]  The  boundary  element  methods  for  optimal  boundary  control  of  elliptic 


partial  differential  equations 

Elliptic  PDEs  often  appear  in  continuum  mechanics.  Assume  that  the 
control  appears  on  the  boundary.  Traditional  numerical  methods  like  the 
finite  differences  or  finite  elements  require  the  discretization  of  the 
entire  domain.  The  boundary  element  method  only  requires  the  discretization 
of  the  boundary  which  yields  numerical  solutions  of  optimal  boundary  control 
much  more  efficiently.  In  [6],  we  treated  the  Neumann  type  boundary  control 
using  the  fundamental  solution.  In  [7],  Dirichlet  boundary  controls  were 
studied  using  the  Poisson  representation  of  harmonic  functions.  Numerical 
solutions  were  computed  which  confirmed  the  theoretical  error  estimates. 


Next,  under  headings  [E]-[G],  we  describe  the  research  that  is  presently 
in  progress. 


[E]  Analysis,  designs  and  behavior  of  dissipative  joints  for  coupled  beams. 


In  [1],  [3],  we  have  studied  the  case  where  a  stabilizer  is  installed  at 
boundary.  What  do  we  know  if  a  stabilizer  is  installed  at  an  in-span  point? 
Let  this  point  be  x  =  a,  0  <  a  <  L.  Assume  that  the  boundary  conditions  at  x 
-  0  and  x  =  L  are  energy  conserving.  Then 


i 


■ji'i  *ji  r.  rr*7*rrrrr>  i  n«  v  m  v  a 
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(15)  ^  [Total  energy  at  time  t]  =  ^  {|  [  my^  ( x .  t )  *  EIyxx  (x ,  t )  ]  dx 

L 

+  \  J  [my^(x.t)  *  EIyxx(x, t ) ]dx) 

=  -[yt(«".t)*EIyxx(oT.t)  -  yt(a".t)  •  E 1 yxx ~ t ) ] 

♦  [yxt  («  •  t )  •  EIyxx(ot~,t)  -  yxt(«"-t)  *  E^xxl**,!)] 

<0,  for  all  t  >  0 . 

We  collect  all  possible  transmission  conditions  below  which  can  cause  energy 
dissipation  in  (15): 

(El)  Continuity  of  bending  moment  and  shear: 

-Ely-xxlc-.t)  =  -E!yxx,«*.t) 

-EIyXXX(<x'-t)  =  -EIyxxx(a*-t) 

discontinuity  of  velocity  (and  displacement)  and  rotation: 

yt  C«  ” .  t )  -  yt(a*.t)  =  •  I-EIyxxx(«-.t)J  *  c} [-Elyxx(«’ . t )  ] 

yxt(«'.t)  -yxt(«'.t)  *  c2["EIyxxx{a't,]  ~  k2[_EIyxx(a  t,) 

(E2)  Continuity  of  displacement  and  bending  moment 

v  ( a  .  t  )  -  y  ( a  *  ,  t ) 

* E 1  y x x ( a  .  t  )  -  -EIyXJt(«*.  t) 
discontinuity  of  shear  and  rotation 

f'EIyxxx(a  •t,l  ‘  I-EIyxxx,*'-t»]  =  *<f  Vt  («  - 1 )  *  Cj  [  -E I  vxx  (a  .  t )  ] 


4 


* 
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yxt(a’’t)  ■  yxt(a  +  lt)  =  c2yt(«'.t)  -  k|[-EIyXJ|(a  .t)]. 

(E3)  Continuity  of  displacement  and  rotation 

y(«", t)  =  y(o+, t) 
yx(«'.t)  =  yx(«+,t) 

discontinuity  of  shear  and  bending  moment 

t~EIyxxx(c,’-t)J  "  f-EIyxxx(o+-t)]  =  -  c1yxt(«“,t) 

-  t-EIyxx(a+ , t) ]  =  c2yt(«-.t)  -  k2yxt (a" , t ) . 

( E4 )  Continuity  of  rotation  and  shear 

yx(«-.t)  =  yx(<*”,t) 

-EIyxxx(a’>t:)  =  ^^xxx^'^) 

discontinuity  of  velocity  (and  displacement)  and  bending  moment 

yt(«“.t)  -  yt(«\t)  =  kj[-EIyxxx(a_,t)  ]  *  ciyxt <«"■*> 

[-EIyxx(«-.t)]  -  t-EIyxx(a".t)]  =  c2t -Eryxxx  (<*"  .  t )  ]  -  k|yxt(a‘,t). 
Note  that  throughout  (E1)-(E4),  k^  >  0,  k2  >  0  and 

kf  *  >  0 

(16)  ( c 1 _c2 *  '  '  k\<2  ~  k2^2  *  °*  for  all  <-V  e  R. 

In  each  of  these  cases,  if  a  physical  variable  is  discontinuous,  then  its 


conjugate  variable  must  be  continuous. 


These  mechanical  designs  included  here  were  worked  out  jointly  with 


Professor  Harry  H.  West  of  the  Civil  Engineering  Department  of  The 
Pennsylvania  State  University,  who  is  specialized  in  structural  analysis  and 
designs . 


Eigenf requency  analysis  along  with  the  aforementioned  designs  will  be 
presented  in  a  forthcoming  paper  [5]. 


[F]  Study  of  a  new  beam  model  proposed  by  D.  L.  Russell 

In  [9],  D.  L.  Russell  proposed  the  following  beam  model 

mwtt(x.t)  *  Elwxxxx(x,t)  -  hi  J^(x.m*xt<x.t)-w(tU.t)]d* 

0  <  x  <  1 ,  t  >  0 


with,  say,  clamped  left  end 


w(O.t)  =  0 
wx(0,t)  =  0 


and  free  right  end 


wxx  ^ 1 ' *  ^  0  j 

-EJw^xd^)  *  /  h(l.f)Iwxt(i.t)  -  wftU.t)]df  =  0. 


The  above  model  seems  to  exhibit  damping  characteristics  of  thin  slender 
beams  more  closely  to  a  realistic  beam. 

So  far,  we  have  completed  the  investigation  of  the  case  when 

h(x,f )  a  constant  =  h. 

Our  preliminary  results  are  the  following, 
i )  When 

n  >  • 


ft 


" overdamping"  occurs. 


ii)  When 

h  _  0  /  m 
E7  '  1 

"critical  damping"  occurs . 

i i i )  When 


"underdamping"  occurs . 


iv)  There  exists  a  holomorphic  Cg-semigroup  corresponding  to  the  dynamic 
beam  model . 

This  research  is  currently  being  continued.  The  results  have  potential 
applications  to  the  modelling  and  control  of  structures,  and  provide  some 
unifying  ground  between  distributed  parameter  and  finite  element  (i.e., 
lumped  parameter)  theories. 

[G]  Numerical  computation  of  the  spectrum  of  a  boundary  damped  dynamic 

blharmonic  equation  modelling  a  plate 

Recently.  J.  Lagnese  [8]  proved  the  exponential  decay  of  energy  of  a 
dynamic  plate  using  certain  boundary  stabilization  scheme. 

A  Ph.D.  student  is  writing  computer  codes  to  program  the  Legendre 
spectral  method  for  the  boundary  damped  plate.  We  expect  to  run  the  program 
at  the  John  von  Neumann  Center  using  Cyber  205.  The  work  involving  vectori- 
zation  and  parallel  algorithms  is  in  good  progress. 


Fig.  5  Stabilizer  arrangement  satisfyin; 
f-EI  yxxx(L,t)  =  -k?yt(L-t)’  kl  >  0  ' 

l"EIyxx(L,t)  -  0. 


kJyt(L,t) 


2 

.  Shear(L.t)  -  kj.  yt(L,t) 
Vertical  Damper 


7  Stabilizer  arrangement 


-EIyxxx(L>t)  =  o 

-EIyxx(L.t)  = 


satisfying 


k|  >  0. 


cd 


yxt(L*t)h/Z 


Moment (L, t) 


k2  >r*t(L,t) 


h2/4 


Horizontal  Damper 


DispliL.t)  -  y(L , t)  -  0 
Moment ( L, c)  -k?y  (Lft)  h2/4 


yxt(L,t)h/2 


Inclined  Damper 


10  Stabilizer  arrangement  satisfyin; 
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